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When a black hole is in an empty space on which there is no matter field except that of the 
Hawking radiation (Hawking field), then the black hole evaporates and the entropy of the black 
hole decreases. The generalised second law guarantees the increase of the total entropy of the 
whole system which consists of the black hole and the Hawking field. That is, the increase of 
the entropy of the Hawking field is faster than the decrease of the black hole entropy. In naive 
sense, one may expect that the entropy increase of the Hawking field is due to the self-interaction 
among the composite particles of the Hawking field, and that the sei/-relaxation of the Hawking 
field results in the entropy increase. Then, when one consider a non-self-interacting matter field 
as the Hawking field, it is obvious that the self-relaxation does not take place, and one may think 
that the total entropy does not increase. However, using nonequilibrium thermodynamics which has 
been developed recently, we find for the non-self-interacting Hawking field that the rate of entropy 
increase of the Hawking field (the entropy emission rate by the black hole) grows faster than the rate 
f> |. of entropy decrease of the black hole along the black hole evaporation in the empty space. The origin 

D ■ of the entropy increase of the Hawking field is the increase of the black hole temperature. Hence an 

understanding of the generalised second law in the context of the nonequilibrium thermodynamics 
is suggested; even if the self-relaxation of the Hawking field does not take place, the temperature 
increase of the black hole during the evaporation process causes the entropy increase of the Hawking 
field to result in the increase of the total entropy. 
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When a black hole is in an empty space on which there is no matter field except that of the Hawking radiation 
(Hawking field) [l| , then an outgoing energy flow from the black hole into the empty space arises due to the Hawking 
radiation and the black hole loses its mass energy. This denotes that a black hole evaporates. Concerning the black 
bi). hole evaporation in the empty space, the generalised second law (GSL) of the black hole thermodynamics is an 
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I. INTRODUCTION 



space, 

important and interesting issue 01 O- The GSL conjectures that the entropy of a black hole S g is given by the 
horizon area (divided by 4) and that the total entropy of the whole system, which consists of the black hole and the 
Hawking field, increases during the black hole evaporation process. Note that, when the GSL holds, the black hole 
entropy S g decreases during the evaporation process, but the entropy of the Hawking field increases faster than the 
- - ■ decrease of S g . 

By the way, the outgoing energy flow from the black hole into the empty space denotes that the Hawking field is in 
a nonequilibrium state, since no energy flow arises in an equilibrium state. Therefore the black hole evaporation in 
the empty space should be interpreted as a nonequilibrium process of the whole system. In order to prove the GSL 
in the context of the black hole evaporation in the empty space, we have to define a nonequilibrium entropy of the 
Hawking field. 

However, because a general definition of the nonequilibrium entropy of matters has not yet been formulated, the 
existing proofs of the GSL consider an equilibrium between a black hole and a heat bath surrounding it, or assume 
the existence of a well-defined nonequilibrium entropy of an arbitrary matter Q. Here note that, in the existing 
proofs, a self-interacting matter field has been considered as the Hawking field. For the black hole evaporation in 
the empty space, the self-interaction causes the se?/-relaxation of the Hawking field itself during propagating on the 
empty space to increase the entropy of the Hawking field (if a nonequilibrium entropy were to be defined). Then one 
may expect naively that the self-relaxation of the Hawking field is the origin of the increase of the total entropy of 
the whole system. This seems to be a very reasonable picture of the GSL in the context of the black hole evaporation 
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in the empty space. However, I dare to suggest an understanding of the GSL in the context of the nonequilibrium 
thermodynamics as follows. The self-interaction of the Hawking field is not necessary, but the increase of the black 
hole temperature during the evaporation process plays an essential role for the validity of the GSL. That is, due to 
the temperature increase, the rate of entropy increase of the Hawking field (the entropy emission rate by the black 
hole) grows faster than the rate of entropy decrease of the black hole, then the total entropy increases without the 
self-interaction of the Hawking field. Note that the increase of the black hole temperature along the loss of mass 
energy is due to the negative heat capacity of the black hole. The negative heat capacity is peculiar to self-gravitating 
systems [j| . Therefore our understanding of the GSL depends strongly on the feature of the self-gravitation of the 
black hole. 

In order to explain our understanding of the GSL, it is appropriate to consider a non-self-interacting matter field as 
the Hawking field. For example, photon, graviton, and a massless free Klein-Gordon field (□<£> = 0) can be the non- 
sclf-interacting Hawking fields. Recently, in the study on nonequilibrium thermodynamics of a non-self-interacting 
matter field, it has been reported that a nonequilibrium thermodynamics for a radiation field (photon gas) was 
formulated in a consistent way @. With referring to report Q, this paper aims to calculate the time evolution of 
the total entropy of a black hole and the non-self-interacting Hawking field in a simplified model of the black hole 
evaporation in the empty space, and to give an understanding of the GSL as explained above. Further, since report 
treats the photons which are collisionless and massless particles, we consider that the Hawking field is not only 
non-self-interacting but also massless. 

This paper is organised as follows. The strategy to reach our understanding of the GSL is explained in section ITT1 
In section HTT1 the nonequilibrium entropy obtained in is briefly explained. Section HVI is devoted to the application 
of the nonequilibrium entropy to our strategy. Summary and discussions are in section IV! 

We set the fundamental constants unity, K = G = c = kg = 1, Then the Stefan-Boltzmann constant is a — 7r 2 /60, 
which is appropriate for a photon gas. When one consider a general non-self-interacting massless field as the Hawking 
field, replace the Stefan-Boltzmann constant as a — > a' — iV7r 2 /120 = N a/2, where N — n\, + (7/8) n/, and rib 
(rib = 2 for photons) is the number of helicities of massless bosonic field and n/ is that of massless fermionic field. 



II. STRATEGY TO REACH OUR UNDERSTANDING OF THE GSL 



A. GSL in the context of the black hole evaporation 



Before proceeding to the explanation of our strategy to reach the understanding of the GSL in the context of 
nonequilibrium thermodynamics, we summarise the statement of the GSL in the context of the black hole evaporation 
in the empty space. 

Some classical properties of a black hole solution of the Einstein equation can be mathematically formulated so 
that they correspond to the four laws of the ordinary thermodynamics On the other hand, the Hawking radiation, 
which is a quantum effect of matter fields on a classical black hole spacetime, implies that the black hole emits a 
thermal radiation of temperature T g = k/2tt, where k is the surface gravity at the black hole horizon Q. Then, 
because of the Hawking radiation, the mathematical correspondence between the properties of a black hole and the 
laws of thermodynamics seems to be physically reasonable. This mathematical correspondence has been known as 
the black hole thermodynamics. In the black hole thermodynamics, the quantity S g which corresponds to an entropy 
of the black hole is determined to be S g = A g /4, where A g is the spatial area of the black hole horizon. 

Hereafter we consider the Schwarzschild black hole for simplicity. Consider the case that the black hole is in an 
empty space on which no matter field except the Hawking field exists, and the black hole evaporates due to the 
Hawking radiation. The black hole thermodynamics suggests that the black hole behaves like a spherical black body 
whose equations of states are 

where R g and E g are respectively the areal radius and the mass energy of the black hole. 

When the Hawking radiation (the quantum effect of matter fields) is not considered, the mass energy E g never 
decrease and the quantity S g never decrease as well. Then, from classical viewpoint, it seems that S g corresponds 
to the entropy of the black hole. However it is obvious that, when the Hawking radiation is taken into account, S g 
decreases during the black hole evaporation process due to the energy loss by the Hawking radiation. Then the GSL 
conjectures that the quantity S g in equations is the true entropy of the black hole, and the total entropy of the 
whole system increases during the black hole evaporation process in the empty space |2|. 



3 



Here we give some supplementary comments on the energetic and entropic properties of the black hole. For the 
first on the energetics, the equations of states JTJ give the negative heat capacity C g of the black hole, 



dEg 

r 9 



This denotes that, when the energy E g is lost, the temperature T g increases. The negative heat capacity is a peculiar 
energetic property to the self-gravitating systems p| . Therefore it is obvious that the equations of states QJ contain 
the self-gravitational effects of the black hole on its own thermodynamic states. Next for the black hole entropy, it has 
already been revealed that, if there is no black hole horizon (R g = 0) on the spacetime, the quantity which corresponds 
to the entropy of the spacetime itself is zero Q [8j]. That is, the quantity which corresponds to the entropy of the 
whole gravitational field is given by the quantity S g in equations (Q. Hence, the energetic and entropic properties of 
a black hole are well encoded in the equations of states QJ. 



B. Strategy to reach our understanding of the GSL 

We have to define the total entropy of the whole system. It has already been shown in reference @ that, for the 
equilibrium state of the black hole with a general self-interacting Hawking field, the total entropy is given by a simple 
sum of the black hole entropy S g with the equilibrium entropy of the general Hawking field. Therefore, in the case of 
the black hole evaporation in the empty space, we assume that the equilibrium case is extended simply so that the 
total entropy Stot of the whole system which consists of the black hole and the nonequilibrium Hawking field is given 
by the simple sum, 

Stot = Sg + S m , (3) 

where S g is given by equations Q and S m is the nonequilibrium entropy of the general self-interacting Hawking field. 
The GSL conjectures the inequality dStot > during the black hole evaporation process. Here note that, because the 
black hole evaporation proceeds very slowly (as seen later in section flV A|l . we approximate the black hole itself is in 
an equilibrium state at each moment during the evaporation process, while the Hawking field is in a nonequilibrium 
state. Hence we use the equilibrium quantity S g as the black hole entropy in equation (JSJ. More discussions on the 
validity of the additivity and of the existence of a well-defined nonequilibrium entropy S m of an arbitrary matter 
are given later in section Ivl 

Because the Hawking field is in a nonequilibrium state due to the energy flow by the Hawking radiation, the problem 
is how to deal with the nonequilibrium entropy S m . Note that, it is appropriate to assume that the nonequilibrium 
entropy of any self-interacting matter increases during a relaxation process, since it is the self-interaction that causes 
the increase of the matter entropy in the context of the ordinary equilibrium thermodynamics. Therefore, if we find 
the increase of the total entropy of the whole system which consists of the black hole and the non-self-interacting 
Hawking field, then it implies that the GSL holds with an arbitrary self-interacting Hawking field. Here one may 
show an example as an objection that, when an inter-stellar gas collapses to form a star, although the net entropy will 
increase, the self-gravitational effect decreases the entropy of the collapsing gas. However as discussed in appendix lAl 
we can expect that this objection is not true of the Hawking field. 

Next turn our attention from self-interaction to the effect by the eaiernaZ-gravitational field on the Hawking ra- 
diation. It should be noted that, even when the self-interaction (including self-gravitation) of the Hawking field is 
ignored, there remain some gravitational interactions between the Hawking field and the black hole, which are for 
example the gravitational red shift, curvature scattering, lens effect and so on. The external-gravitational field works 
as if a virtual medium on which the Hawking field propagates, then the composite particles of the Hawking field 
interact with the virtual medium to result in the relaxation of the Hawking field toward some equilibrium state. 
Therefore, it is appropriate to assume that the external-gravitational interaction causes the increase of the entropy 
of the Hawking field. Hence we expect that, if we can calculate a nonequilibrium entropy of the non-self-interacting 
Hawking field by ignoring the external-gravitational interaction, then the rate of increase of such an approximate 
entropy is smaller than that of the nonequilibrium entropy of the non-self-interacting Hawking field calculated by 
including the external-gravitational interaction. 

From the above discussions, we can expect that the following inequality holds along the black hole evaporation 
process in the empty space, 

dS^t > dS\T ] , (4) 

where Stot is the total entropy with a general self-interacting Hawking field given by equation which includes 
the gravitational interactions between the black hole and the Hawking field, and S^™"^ is the total entropy of the 



4 



whole system which consists of the black hole and the non-self-interacting Hawking field. Here the entropy S t ™ t m is 
decomposed as 

Q (min) _ c , q (z\ 

where S ns i is the nonequilibrium entropy of the non-self-interacting Hawking field without including the gravitational 
interaction with the black hole. Therefore, if the inequality dSj%f n ^ > holds, then we can expect that the GSL, 
dStot > 0, follows. 

In naive sense, one may think that the nonequilibrium entropy S ns i of the non-self-interacting Hawking field does 
not increase, because the self-interactions and the external-gravitational potential are ignored and the Hawking field 
never relaxes by itself. The absence of the self-relaxation is a right conclusion from above discussions. However it 
should be emphasised that there still remains an interaction, the emission process of the Hawking radiation by the 
black hole. Here note that the negative heat capacity C g of the black hole denotes the increase of the black hole 
temperature T g along the loss of the mass energy E g . Therefore it can be expected that the emission process causes 
the increase of the entropy Snsi of the non-self-interacting Hawking field even when the gravitational interaction with 
the black hole is ignored. The higher the temperature T g grows, the more amount of matter entropy is emitted with 
the Hawking radiation. 

In this paper, by constructing a concrete formula of the nonequilibrium entropy S ns i of a non-self-interacting 
Hawking field without including the gravitational interaction with the black hole, we aim to show that the inequality 

dSt™ n ^ > holds during the black hole evaporation in the empty space. The concrete method for the calculation of 
the entropy S ns i is explained in the following sections. 



C. The model of the black hole evaporation 

Because we are interested in the non-self-interacting Hawking field without the gravitational interaction with the 
black hole, it is enough to consider the model that a black body satisfying the equations of states Q is put in a flat 
spacetime. This modelling retains the self-gravitational effects of the black hole on its own thermodynamic states 
(negative heat capacity), as mentioned in section III Al Therefore, in order to pick up the nonequilibrium effects 
of the Hawking field on the black hole evaporation process, we consider a simplified model named NHR after the 
nonequilibrium Hawking radiation: 

NHR Model: Consider a spherical black body satisfying equations (|TJ, which emits some non-self- interacting mass- 
less matter fields with Planckian spectrum (free Hawking fields). Then put the black body in a flat spacetime 
on which there exists no matter field except the free Hawking fields. Then the black body evaporates due to 
equations and the emission of the free Hawking fields. This evaporation can be understood as a relaxation 
process of the whole system which consists of the black body and the fee Hawking fields. 

Here note that, although the self-interactions (including the self-gravitation) of the Hawking field are ignored, however 
the self-gravitational effects of the black hole on its own thermodynamic state is encoded in the equations of states 
l[T|l . Henceforth we call the black body the black hole, and call the non-self- interacting massless fields the free Hawking 
fields. 

In the NHR model, the black hole continues to emit the free Hawking fields without any in-falling energy into the 
black hole, and the free Hawking fields continue to spread out into the empty space. That is, an outgoing energy flow 
from the black hole into the empty space exists during the black hole evaporation process. The outgoing energy flow 
equals the rate of energy increase in the free Hawking fields. 

Further because of the quasi-equilibrium approximation of the black hole itself (see the first paragraph in section 
111 Hp . wc notice that, from the thermodynamic viewpoint, the thermodynamic state of the black hole itself during 
the evaporation process changes along a sequence of equilibrium states characterised by equations Q in the state 
space, and that the black hole behaves always as a perfect black body of temperature T g at each moment during the 
evaporation process. (Recall the quasi-static process in the ordinary equilibrium thermodynamics.) Therefore the 
rate of energy increase J in the free Hawking fields in the NHR model is given by the Stefan-Boltzmann law, 

J = a'T g i A g , (6) 

where a' is the general Stefan-Boltzmann constant given at the end of section H] 

This J is equivalent to the energy flow in the free Hawking fields, and denotes that, while the thermodynamic state 
of the black hole changes along a sequence of equilibrium states, the thermodynamic state of the free Hawking fields 
should change along a sequence of nonequilibrium states in the state space, where each nonequilibrium state in the 
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FIG. 1: MS model. The radiation field between the outer and inner black bodies is in a steady state. The entropy of the 
radiation field is given by the steady state thermodynamics for a radiation field. 

sequence of the field's states should have the energy flow J. Hence the time evolution of the total entropy 
should be analysed with a nonequilibrium thermodynamics for the free Hawking fields. In order to carry out the 
nonequilibrium analysis, we refer report [|| which has formulated a nonequilibrium thermodynamics for a radiation 
field (photon gas) in a consistent way. This can easily be extended to a general non-self-interacting massless field as 
shown in the next section IIIII Then the nonequilibrium entropy S ns i of the free Hawking fields is concretely given 
using the formula obtained in |fj. 

Here one may point out that, for an equilibrium of a black hole with the Hawking field, the gravitational red shift 
results in the need for a higher temperature of the Hawking field than the black hole temperature T g in order to retain 
the equilibrium against the attractive gravitational force by the black hole . Then an energy flow stronger than J 
would be appropriate in order to retain the effect of the gravitational red shift. However the NHR model ignores the 
external-gravitational potential, and the energy flow J of equation © is suitable for our purpose. Further note that a 
lower temperature of the free Hawking fields implies a smaller emission rate of entropy by the black hole. Therefore, 
by ignoring the gravitational red shift, we can recognise that inequality J3J becomes more reasonable. 

Although some simplifications and assumptions are put in the NHR model, which one may think as a too simplified 
model, we will use it to obtain some insight into the GSL. Further in order to avoid confusions, it is helpful for readers 
to mention here that the so-called black hole phase transition (or the Hawking-Page phase transition for AdS black 
holes) in the ordinary black hole thermodynamics 01 is ou t of the scope of this paper. The black hole phase transition 
is the conclusion which comes from the study on the equilibrium states of the system in which a black hole is put in 
a heat bath. Here in this paper we do not consider any heat bath surrounding the black hole, but put a black hole in 
an empty and infinitely large space where the black hole can not be in any equilibrium with its environment. We are 
interested in the dynamical process of a black hole evaporation in the empty space as a nonequilibrium process, and 
aim to clarify the effect of the increase of black hole temperature on the GSL. 

III. STEADY STATE THERMODYNAMICS FOR FREE HAWKING FIELDS 

As a middle step before proceeding to the calculation of the nonequilibrium entropy S ns i of the free Hawking fields 
in the NHR model, we consider a more simplified model which considers a radiation field (photon gas) as the free 
Hawking fields and can be realised in a laboratory experiment. This simple model is shown in figure ^ an d named 
MS after the middle step: 

MS model: Make a vacuum cavity in a large black body of equilibrium temperature T out and put an another smaller 
black body of temperature Tj„ (^ T out ) in the cavity. For the case Tj„ > T out , the radiation field (photon gas) 
emitted by the two black bodies causes an energy flow from the inner black body to the outer one. When the 
outer black body is isolated from the outside world, the whole system which consists of the two black bodies 
and the radiation field between them relaxes to a total equilibrium state in which the two black bodies and the 
radiation field have the same equilibrium temperature. 

Note that the equations of states of the two black bodies are not specified. Here following report [||, we consider that 
the black bodies emit a radiation field (photon gas). However it is easily extended to a general non-self- interacting 
massless field by replacing the Stefan-Boltzmann constant a with the generalized one a 1 as mentioned at the end of 
section^] In section llVl we will extend the MS model to the NHR model. To do so, we will make the inner black body 
have the equations of states and remove the outer black body by setting the cavity infinitely large and T out = 0. 

Concerning the relaxation process required in the MS model, consider the case satisfying the following two condi- 
tions. The first is that the thermodynamic state of each black body evolves so slowly that it changes along a sequence 
of equilibrium states in the state space during the relaxation process. The second condition is that the volume of 
the cavity is so small that the speed of light is approximated as infinity during the relaxation process. These two 
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conditions denote that [6(, at each moment of the relaxation process, the retarded effect of photons is ignored, and the 
thermodynamic state of the radiation field sandwiched by the two black bodies is well approximated by the nonequi- 
librium state which possesses a stationary energy flow J ra d determined by the temperatures of the two black bodies 
Ti 7l and T out . The energy flow J rao i is given by the Stefan-Boltzmann law, J ra d = a ( T^\ — T 4 ut ) Aj„ , where Ai n is 
the surface area of the inner black body. This J ra d equals the energy exchange par a unit time between the two black 
bodies via the radiation field. 

From above discussions, it is enough for us to consider a set (state space) of the steady states (macroscopically 
stationary nonequilibrium state) of the radiation field whose nonequilibrium nature is raised by the steady energy 
flow J r ad- That is, the sequence of the nonequilibrium states of the radiation field during the relaxation process lies in 
the set (state space) of the steady states. Although each steady state alone does not cause any time evolution of the 
system, however the time evolution of the temperatures T in and T out can be obtained by shifting the thermodynamic 
state of the radiation field along the sequence of the steady states as time goes on. Therefore we need a consistent 
definition of the steady state entropy of the radiation field. Once the steady state entropy of the radiation field is 
obtained, we can express explicitly the total entropy during the relaxation process of the whole system which consists 
of the two black bodies and the radiation field. Further, by extending the definition of the steady state entropy into 
the case of an infinitely large cavity with T out = 0, it is expected that the total entropy of the NHR model can 

be analysed explicitly as well. 

A consistent thermodynamic framework for the steady states for a radiation field in the MS model has already been 
constructed 0- The outline of the construction of the entropy in the steady state thermodynamics for a radiation 
field in the MS model is as follows. 

Refer the Landau-Lifshitz type definition of a nonequilibrium entropy for a bosonic gas pl| , 

f dp 3 

Sboson = J 7^)3 9 ^ 3 ^ 1 + ) 111 ( 1 + fa' 3 ) ~~ fP' S 111 ] ' ( 7 ) 

where p is the momentum of a photon, x is a spatial point, and gpg and fpg are respectively the number of states 
and the average number of particles at a point (p, x) in the phase space of the bosonic gas. Note that it has also been 
shown in reference that the maximisation of Sboson for an isolated system (SSboson — 0) gives the equilibrium 
Bose distribution. This is frequently referred to in many works on nonequilibrium systems as the H-theorem. However 
in reference concrete forms of gpg and fpg are not specified, since an arbitrary system is considered. 

In the MS model, the radiation field is sandwiched between two black bodies. Then we can determine gg t g and fg$ 

as 

gps — 2 (helicities of a radiation field) , fgs = f — ,„ , _ =^ , (8) 

exp[o;/J (p 7 x)\ — 1 

where the frequency of a photon uj — \p\, and T{p, x) is given by 

mt - ->n J T in for p = pin at X 

1 (p, x) = < , 
\T out for p = pout at x 

where pi n is the momentum of the photon emitted by the inner black body and p ou t by the outer black body. The 
x-dependence of T(p, x) arises, because the directions in which the photons of pi n and p ou t can come to a point x 
vary from point to point. 

From above, we obtain the steady state entropy of a radiation field, S ra di 

Srad = J dx 3 s ra d(x) , s rad (x) = ( g in (x) T in + g ou t (x ) T Q \ t ) , (9) 

where gi n (x) is the solid angle (divided by 47r) covered by the direction of Pi n at x as shown in figure ^ an d 9out{x) 
is similarly defined with p ou t- In deriving JSJ, we used the relations, J °° dx [x 2 /(e x — 1)] ln(e x — 1) = ll7r 4 /180 and 
/ °° dx [x 2 /{l - e- x )] ln(l - e-*) = -^ 4 /36 . 

Further, with defining the steady state internal energy E ra d = 2 J dp 3 dx 3 uj fp t $ and the other variables such as the 
free energy with somewhat careful discussions, it has already been checked that the 0th, 1st, 2nd and 3rd laws of 
the ordinary equilibrium thermodynamics are extended to include the steady states of a radiation field. That is, the 
steady state thermodynamics for a radiation field has already been constructed in a consistent way. Especially in the 
steady state entropy 10, it has already been revealed in [(| that the total entropy of the whole system which consists 
of the two black bodies and the radiation filed increases monotonically during the relaxation process required in the 
MS model, dSi n + dS ou t + dS ra d > 0, where the equality holds for the total equilibrium state, Si n and S ou t are the 
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entropies of the two black bodies which are given by the ordinary equilibrium thermodynamics, and it is assumed 
that the two black bodies are made of ordinary materials and have positive heat capacities. 

When we consider a massless fermionic field instead of the radiation field in the MS model, equation (JJJ) is replaced 

by EH 

f dp 3 3 

&f ermion — / To vT^^ 9p,x [ ( 1 fp,x ) In ( 1 fp,x ) fp.x hi ffj^x ] > (10) 
J (Z-K)- 3 

where g^g = rif is the effective number of helicities, and f^g = [exp[w/T(p, x)] + 1] — 1 . Then, following the same 
procedure as a radiation field, we can obtain the steady state entropy of a massless fermionic field as 

Sf = 8 ~~2 Srad ' ^ 

where we used, J™ dx [x 2 /{e x + 1)] lmV + 1) =tt 4 /60 + A, and / °° dx [x 2 /{l + e- x )} \n{l + er x ) = 11tt 4 /180-A, 
where A = [(ln2) 2 - tt 2 ] (ln2) 2 /6 + 40(4, 1/2) + (7 In 2/2) C(3), CO) is the zeta function and 4>{z,s) is the modified 
zeta function (Apell's function). Hence, when the black bodies in the MS model emit rif massless fermionic modes 
and rib massless bosonic modes (rib = 2 for a radiation field), the steady state entropy of the matter fields between 
the two black bodies is given by equation by replacing a by a' given at the end of section [I] 



IV. TIME EVOLUTION OF THE TOTAL ENTROPY 



A. From the MS model to the NHR model 

We proceed to the analysis of the total entropy >S{™ in the NHR model, which can be obtained from the MS 
model by considering the equations of states of the inner black body as equations (0 and by removing the outer black 
body. The way of removing the outer black body from the MS model is to set the cavity's volume infinite and to 
set the outer temperature zero T out = 0. In the NHR model, the free Hawking fields continue to spread out into the 
empty space. 

The time evolution of the black hole radius R g is given by equations and © , 

— =-J => RM) = R (l ^-_V /3 = L_, (12) 

where Rq = R g {0) is the initial radius, and it is assumed that the emission of the Hawking radiation starts at t = 0. 
This R g {t) leads the time evolution of the thermodynamic state of the black hole and the steady state of the free 
Hawking fields. Equation i|12|) gives the evaporation time (life time) of the black hole, 

R 3 

U = 1280^. (13) 

On the other hand, the time scale t^h that a massless particle travels across the size of the black hole is given by 
tbh = Rq- This time tbh gives a typical time scale of the free Hawking fields to spread out into the empty space. 
Therefore, when the inequality t ev /tbh ^ 0(1) holds, the black hole evaporation proceeds very slowly in comparison 
with the time scale tbh- Henceforth, we assume the order relation t ev /tbh » 0(1), which denotes that the quasi- 
cquilibrium approximation of the black hole itself is valid as mentioned in section III Bl More discussion on this 
assumption is given later in section Ivl 



B. Nonequilibrium entropy of the free Hawking fields 



In order to apply the steady state thermodynamics [a| to the free Hawking fields in the NHR model, the distribution 
function of photons (jHJ should be modified to include the retarded effect of composite particles (rib bosons and rif 
fermions) travelling in an infinitely large space. In this section, we construct the modified distribution function for the 
first, then obtain the nonequilibrium entropy of the free Hawking fields in the NHR model by substituting the modified 
distribution function into formulae 10 and llOjl . Here note that, because of the quasi-equilibrium approximation of 
the black hole itself, we ignore the special relativistic Doppler effects due to the shrinkage of the surface of the black 
hole. 
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R g (t)-~-' 



(t , r , 9 P ) 



FIG. 2: Retarded effect on the Hawking radiation. 



Hereafter, r denotes the areal radius from the centre of the black hole. In order to take the retarded effect into 
account, consider a massless particle which has been emitted at the surface of the black hole at time t and reaches 
the spatial point of r at time t (> t). The emission time t depends not only on the coordinates (t, r), but also on 
the angle 9 p between the radial direction of r and the momentum of the particle (see figure [2J. This emission time, 
t(t, r, 9 P ), can be obtained as a root of the following equation, 



R g {tf = {t-i) + r 2 - 2 (t- t) r cos9 p . 



(14) 



Note that this is the equation of degree six about t, and the appropriate root as the emission time t(t, r, 9 P ) is the 
maximum root in the range, < t < t. Although another root may exist in this range, however such a root corresponds 
to the particle emitted at point q in figure [21 which is obviously unphysical. And the other four roots of this equation 
may be complex valued ones. Further it is also obvious that the angle 9 P has an upper bound, < 9 p < 9 m (t, r). This 
upper bound is given by appendix IB! as (see equations (|B1|) . 1B2|) and l|B3|) ) 



cos 9 m (t,r) 



cos 9 m {t, r) = 



t - t m (t, r) 



t 2 + r 2 



2tr 



for r <R 

, for r > i? and t < r — R 

, for r > Rq and r — R < t < ^/r 



R 



- [t - t m (t, r) } , for r > R and y/ V 2 - R 2 , < t 



where t m (t,r) is the real valued root in the range, < t rn < i, of the equation, 

r 2 - R g {t m f + (t-t m ) 2 . 
Then the distribution function of the free Hawking fields is given as 



1 



f(t,r;u,O p ) 



exp 




io/f 



± 1 



(15) 



for 9 p <9 m (t,r) 
for 9 P > 9 m (t,r) 

where T = T g (t(t,r,9 p ) ), co = \p\, and the signatures "— " and "+" are respectively for bosons and fermions. Here 
note that, while the t- and r-dependence of / expresses the spacetime dependence, the to- and # p -dependence denotes 
the dependence on the momentum p of particles of the free Hawking fields. 

The entropy, S ns i(t), of the free Hawking fields can be obtained with substituting the distribution function l|15|) 
into the place of fp^ in formulae J7J) and (|10|l . 



t+Ro 



drAnr 2 s nsi 



(t,r) 



N 



2880 7T y Vm(t)P) uv R g (t(t,r,y p )) 3 ' 



(16) 



where the same integrals used in obtaining equations 10 and (|llfl are also used, and we set y p = cos# p , y m — cos9 m 
and gpg — nj, and nj for bosons and fermions respectively. Here note that, since we assume that the emission of the 
free Hawking fields starts at t = 0, the free Hawking fields fill the space in the range, R g (t) < r < t + Rq . 
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FIG. 3: Time evolution of the normalised total entropy £j™"'(t) = S 1 ^™ (t)/^™ , where r is the time normalised by the 
evaporation time t ev . The plotted curves are converging as A (= t ev /Ro) increases. 

It is useful to consider here the validity of the nonequilibrium entropy (|16|l . We pay attention to the quasi-equilibrium 
approximation of the black hole itself. This approximation require, not only that the relation t ev /tbh 3* 0(1) holds, 
but also that the shrinkage speed of the surface of the black hole is less than unity, 



dR g (t) 



dt 



-fi 



3 t eV Rg 



< 1 



t < t e 



N 



48\/l57r 



(17) 



This denotes that, when N < 48 2 x 157T ~ 108573, our entropy (|16|) is valid at least until one Planck time before the 
evaporation time t ev . 



C. Total entropy 

From the above discussions, the total entropy of the NHR model which consists of the black hole and the free 
Hawking fields is expressed as S^ m ^ (t) = S g (t) + S ns i (t) , where S g (t) is given by equations (Q and , and S ns i (t) 
by equation (|16|) . In order to show the numerical plot of this total entropy, we normalise as follows, 

T-— T-— T-— I W-M-fl-r^ y (min) , s _ S[™ t m) (t) _ S nsi (t,r) 

tev Tev tt Ho £V (Q) ^ 0) 



Then the normalised total entropy is 



16 / , / , x 



slrV) = X g (rf + ^ / dx dy p ——f- — , (18) 

9A7 Xg(r) Jy m (r.x) X g {T{T,x,y p )Y 



where A is given by 



A = 1280 7T — = ^ = ^ 
N R tbh 



Here note that, because the entropy E^*^ does not explicitly depend on i?o (or N), then, once the value of A is 
fixed, the choice of i?o (or N) is arbitrary with adjusting the value of N (or Rq) appropriately to match with A. And 
we assume A 3> 0(1) as mentioned in section ITY Al 

The numerical results are shown in figures |3] and 0] The parameters i?o and N are chosen so that the relations 
A ^> 0(1) and N < 48 2 x 157T are satisfied, and that the time interval < r < 0.999 is included in the range i|17[l . 
The numerical plots are carried out in the interval < r < 0.999 with the Mathematica version 5.2. 
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FIG. 4: Spatial and temporal evolution of the entropy density a nS i(T,x) and its spherical shell density 4ttx 2 a n si(r, x) for 
(A —)t ev /Ro = 10 4 , where a nS i(T,x) = s nS i (t, r)/£^™ n ' (0) is the normalised spatial density of the steady state entropy, 
x — r/Ro is the normalised spatial distance and r = t/t ev is the normalised time. The arrows denote the direction of null world 
lines (geodesies). 

Figure shows the time evolution of the total entropy £^ n) for A = 10,50, 100, 10 3 , 10 4 , 10 5 and 10 6 H3- The 
plotted curves are converging as A increases. The plots for A > 10 3 are almost coincident. What we can find with this 
figure is as follows. 

• Because the plotted curves are converging, it is appropriate to conclude that the total entropy is monotone 
increasing for A 3> 0(1), and the GSL is well supported. 

• Since the normalised black hole entropy X g (r) 2 is obviously monotone decreasing, the entropy of the free 
Hawking fields increases monotonically faster than the decrease of the black hole entropy. 

• It is suggested rather well that £(™ n ' (r) gives almost the same curves in £-r graph for all values of sufficiently 
large A > 10 3 , and the final value is approximately given by £^"^(0.999) ~ 1.33. 

The left panel in figure 21 shows the normalised spatial density of the entropy of the free Hawking fields, cr ns i{r, x), 
and the right panel shows its normalised spherical shell density, 4irx 2 a ns i(T,x). Figure 0] is made with A = 10 4 , 
however the same behaviour is obtained for all values of sufficiently large A > 10 3 . The arrows in figure 01 denote the 
direction of the null world lines of massless particles. What we can find with this figure is as follows. 

• While the entropy of the free Hawking fields spreads out to decrease the spatial density (left panel) , the spherical 
shell density on each spherical shell remains constant while it spreads out into the empty space (right panel). 
This is consistent with the absence of the self-relaxation in the free Hawking fields. 

• The entropy of the free Hawking fields at the surface of the black hole, which equals the amount of matter 
entropy emitted by the black hole at each moment, increases as the evaporation process proceeds. Here note 
that the the black hole temperature increases obviously. That is, the understanding of the GSL suggested in 
section [J can be extracted from figures and 0] 

V. SUMMARY AND DISCUSSIONS 

When a black hole evaporates due to the Hawking radiation, the thermodynamic state of the Hawking field is in a 
nonequilibrium state. Then, referring to a recent development on the steady state thermodynamics for a radiation field, 
we have analysed explicitly the time evolution of a black hole evaporation process within the framework of the NHR 
model, and obtained the conclusion that the total entropy S^ m ^ increases monotonically for (A =)t ev /Ro 3> 0(1), 
where t ev is the evaporation time and Rq is the initial black hole radius. Therefore, as discussed in section^ we can 
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obtain the understanding of the GSL in the context of the steady state thermodynamics as follows. The increase of 
the black hole temperature due to its negative heat capacity gives the increase of the entropy of the Hawking fields, 
and the rate of entropy increase of the Hawking fields is faster than the rate of entropy decrease of the black hole. 
That is, the self-interaction of the Hawking fields is not necessary for the validity of the GSL. 

Here we pay attention to the last stage of the evaporation process in the NHR model; the total entropy Stot 
increases in an accelerating fashion as shown in figure According to the relaxationprocess in the MS model which 
reaches a total equilibrium state as mentioned in section IIIII or shown in report |(| in detail, one may think the 

accelerated increase of is strange, because the black hole evaporation process in the framework of the NHR 

model is also described as the relaxation process of the whole system. However the accelerated increase of £ t ™ 
never contradict the discussions in section ITTll and in report @, because the black hole has the negative heat capacity 
as shown by equation (J2)l . The negative heat capacity is a peculiar property to the self-gravitating systems Q ■ While 
the positive heat capacity of an ordinary system makes the entropy emission rate decrease to zero as the relaxation 
process proceeds, however the negative heat capacity of the black hole makes the entropy emission rate increase as 
the black hole evaporation process proceeds as shown by the right panel in figure 0| This increase of the entropy 
emission rate is the essence of our understanding of the GSL. 

Here it is recalled that the negative heat capacity of a self-gravitating gas system (not of a black hole) causes 
the accelerated increase of its temperature, which is called the gravothermal catastrophe |5|. The gravothermal 
catastrophe is the phenomenon suggested by the ordinary equilibrium thermodynamics and statistical mechanics. 
Hence the accelerated increase of the total entropy of the black hole and the Hawking fields shown in figure may be 
interpreted as a nonequilibrium version of the gravothermal catastrophe. 

In this paper, the monotone increasing nature of the total entropy is concluded by the numerical results 

shown in figure 01 Therefore, correctly speaking, our conclusion is not analytically completely supported at present. 
If the inequality dSj:™ / dt > will be obtained analytically, our understanding of the GSL will have more rigid 
mathematical reasoning. Further our discussion given in this paper may be understood as a new proof of the GSL, 
dStot > 0, according to inequality (@J. 

We turn our discussion to the validity of the assumption A = t ev /tbh 3> O(l), which guarantees the validity of the 
quasi-equilibrium approximation of the black hole. In semi-classical treatment, the black hole temperature is less than 
the Planck temperature. Then the exotic particles are not expected to be created in the Hawking radiation, and only 
the ordinary well-known standard particles will be created. Hence we can expect that the number of species of the 
Hawking fields (~ N) is less than 108573 = O(10 5 ) (see inequality (|T7J> ). For the case N ~ O(10 5 ), the assumption 
A ^> 0(1) is valid for a black hole of initial radius Rq > O(10) in Planck unit. And if we estimate the order of N 
by the ordinary standard particles (three generations of quarks and leptons and their inner states, gauge particles of 
four fundamental interactions and their helicities), then we find N ~ O(10) and the assumption A 3> O(l) is valid for 
a black hole of initial radius Rq > O(10 _1 ) in Planck unit. Hence, when no exotic particle appears in the Hawking 
radiation, our discussion based on the NHR model seems to be valid for a semi-classical black hole of initial radius 
Rq > 1 until one Planck time before the evaporation time t ev . 

Here we notice that, for the case N ~ 0(10), the order relation A > O(10 3 ) corresponds to the black hole of 
initial radius R > 0(10). Therefore we find from figure |3| that, for the black hole of R > 0(10) A > O(10 3 )), 

the normalised total entropy (t) evolves in almost the same increasing fashion and reaches almost the same 

final value ~ 1.33. Further it should be noted here that, even if we consider the very last stage of the evaporation 
process where the quasi-equilibrium approximation of the black hole breaks down and the evaporation may proceed 
to a quantum gravitational process, it is reasonable to expect that the total entropy after the quantum evaporation 
process finished should be greater than the total entropy before the evaporation process proceeds to the quantum one, 
since the well-defined entropy has to be a non-decreasing quantity. Hence, because inequality (0J indicates that S\^ m ^ 
is the lower bound of the total entropy of the whole system, it is suggested that, when the black hole evaporates in 
the empty space with emitting any self-interacting Hawking field, the final value of the total entropy Stot (not the 

lower bound S^ in ^) should be larger than 1.33 x S g (0). 

Finally we discuss about the strategy to reach the understanding of the GSL in the context of the steady state 
thermodynamics. The strategy explained in section^has two key assumptions remained to be proven, the additivity 
of the entropy @ and the existence of a well-defined nonequilibrium entropy S m of an arbitrary self-interacting 
Hawking field. However, these two key assumptions seem to be reasonable as are discussed in the following two 
paragraphs. 

On the validity of the additivity J3j, we refer the present status of the study on the laboratory system which has 
self-interactions but not self- and external-gravitational interactions. For example, consider that the system is in a 
nonequilibrium state which is not so far from an equilibrium state that the heat flux in that system is not extremely 
strong. Then the so-called extended irreversible thermodynamics [T^ | gives a well-defined nonequilibrium entropy up 
to the second order in the expansion by the heat flux. Further the nonequilibrium entropy satisfies the additivity as in 
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the ordinary equilibrium thermodynamics. Therefore, we can expect that, when the equilibrium entropy satisfies the 
additivity even with the gravitational effects, the nonequilibrium entropy also satisfies the additivity. Hence, due to 
the additivity of the equilibrium entropies of the black hole and the Hawking field , it seems reasonable to assume 
the additivity J3J). 

Next we consider the validity of the existence of a well-defined nonequilibrium entropy S m of an arbitrary self- 
interacting Hawking field. As mentioned in the previous paragraph, the nonequilibrium entropy has already been 
defined well up to the second order in the expansion by the heat flux 12]. Further, on the other hand, the evaporation 
time t ev of a black hole has an extremely long time scale, A = t ev /tbh 3> 0(1), as discussed above. Therefore it seems 
appropriate that the black hole evaporation proceeds so slowly that the nonequilibrium state of the self-interacting 
Hawking field is not so far from an equilibrium state that the expansion of the thermodynamic quantities by the 
heat flux up to the second order can be a good approximation. Hence it is reasonable to assume the existence of a 
well-defined nonequilibrium entropy of an arbitrary self-interacting Hawking field. 



APPENDIX A: BLACK HOLE EVAPORATION AND STAR FORMATION PROCESS 



This appendix discusses in a simple way about the similarity between the black hole evaporation process and the 
star formation process. When an interstellar gas collapses to form a star, it is commonly believed that there arises 
the increase of the net entropy of the system which consists of the collapsing gas and the radiated matters from the 
gas. However the se//-gravitational effect of the collapsing gas causes the decrease of the entropy of the collapsing gas 
itself. It is briefly explained as follows. 

Since the pressure of the gas at its surface is zero, the loss of energy AE (par a unit time) of the collapsing gas by 
the radiation of energy is actually the loss of heat due to the first law of the thermodynamics, 

AE ~ T AS (Al) 

where T is the temperature of the collapsing gas, and AS is the loss of entropy (par a unit time). The energy loss 
AE can be interpreted as the minus of the luminosity L of the collapsing gas, L — —AE. Then, with the assumption 
of the local mechanical and thermal equilibrium of the collapsing gas at each moment of its collapse, relation (|A1|) 
can be rewritten more exactly to the following form 

dt 

where U is the total internal energy of the collapsing gas, and O is the total self- gravitational potential given by 

il = dm 



G m 



r(m) 



where G is Newton's constant, M is the mass of the collapsing gas, and the radial distance from the centre of the gas 
r(m) is expressed as the function of mass m inside the sphere of radius r. 

From above, it is recognised that the radius r(m) becomes smaller as the gas collapses, then the self-gravitational 
potential SI decreases to result in the radiation of energy L > and the loss of the entropy AS < 0. Therefore, if 
we consider not a collapsing but an expanding self-gravitating gas, the radius r(m) increases to result in the increase 
of entropy AS > 0. Here note that the Hawking field in the black hole evaporation process corresponds not to the 
collapsing gas but to the the expanding matters radiated by the collapsing gas. Hence it seems reasonable to assume 
that the entropy of a self-gravitating Hawking field is larger than that of a non-self-gravitating Hawking field. 



APPENDIX B: UPPER BOUND OF 9 P 

This appendix derives the upper bound 9 m of the angle 9 p which appears in equation l|14|) . In order to find the 
explicit expression of 9 mi look at the photons emitted by the black hole at the initial time t = 0, which we call the 
initial photons. Obviously, we have 6 m = for t < r — Rq, since no photons has been reached the point. The initial 
photons emitted in radial direction form the boundary of the region which the free Hawking fields fill. The initial 
photons emitted in off-radial directions propagate behind the initial photons emitted in radial direc tions. T herefore, 
the initial photons reach a spatial point of radial distance r > Rq in the time interval, r — Ro<t< y/r 2 — Rq. Within 
this time interval, the upper bound 9 m is given by the initial photon, which is obtained with setting t = in figure 
Then we find 

cos 9 m (t, r) = 1 + 2 t ~ R ° . for r > i? and r - i? < t < \J r 2 - R 2 . (Bl) 
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FIG. 5: Upper bound of 6 V for r < Ro, or r > Ro and t > \/r 2 — R%. 



Further for t > \/r 2 — R\, the upper bound 9 m is not given by any initial photon but by the photon emitted at the 
point b at the time t m shown in figure [S] Then we find 

cos m {t, r) = t JZln^l ; for r > R and yV 2 - R 2 < t , (B2) 

where the time t m (t, r) is a real valued root of the equation, r 2 = R g ( t m ) 2 + ( t — t m ) 2 . This is the equation of degree 
six about t m . The appropriate root for t m should be in the range, < t m < t. This root may be degenerated, since 
the trajectory of this photon is tangent to the sphere of radius R g (t m ), and the other four roots may be complex 
valued ones. Finally turn to the spatial point of radial distance r < Ro. It is obvious for this spatial point that the 
upper bound m is also given by formula l|B2|l. 

cos 6 m (t,r) = t for r < R . (B3) 
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